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O ■ ABSTRACT 

We report the results of a local stability analysis for a magnetized, gravitation- 
ally stratified plasma containing cosmic rays. We account for cosmic-ray diffu- 
^ I sion and thermal conduction parallel to the magnetic field and allow /? = Snp/B'^ 

to take any value, where p is the plasma pressure and B is the magnetic field 
^ I strength. We take the gravitational acceleration to be in the — ^-direction and the 

I equilibrium magnetic field to be in the ^/-direction, and we derive the dispersion 

O ' relation for small-amplitude instabilities and waves in the large- 1 /c^; I limit. We use 

00 ' 

Q I the Routh-Hurwitz criterion to show analytically that the necessary and sufficient 

criterion for stability in this limit is nksdT /dz + dp^/dz + {l/87i)dB'^ /dz > 0, 



where T is the temperature, n is the number density of thermal particles, and 



H ' Per is the cosmic-ray pressure. We present approximate analytical solutions for 

the normal modes in the low- and high-diffusivity limits, show that they are 
consistent with the derived stability criterion, and compare them to numerical 
results obtained from the full, unapproximated, dispersion relation. Our results 
extend earlier analyses of buoyancy instabilities in galaxy-cluster plasmas to the 
/? < 1 regime. Our results also extend earlier analyses of the Parker instabil- 
ity to account for anisotropic thermal conduction, and show that the interstellar 
medium is more unstable to the Parker instability than was predicted by previous 
studies in which the thermal plasma was treated as adiabatic. 



Subject headings: convection — cooling flows — galaxies:active — galaxies:clusters:general — 
magnetic fields — turbulence 
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1. Introduction 

Convection plays an important and well-known role in the transport of energy in stellar 
interiors. It has also been argued that convection is important in a number of low-density 
astrophysical plasmas, such as the intracluster medium in clusters of galaxies (Chandran 
& Rasera 2007) and accretion flows onto compact objects (Quataert & Gruzinov 2000). 
Although convection in stellar interiors has been thoroughly studied over the course of sev- 
eral decades, the theory of convection in low-density plasmas is still being developed, and 
investigations carried out during the last several years have led to some interesting surprises. 

For many years, it was widely assumed that the convective stability criterion for a 
low-density, non-rotating, weakly magnetized plasma is the Schwarzchild criterion, ds/dz > 
0, where s is the specific entropy of the plasma and the gravitational acceleration is in 
the —z direction. However, Balbus (2000, 2001) showed that even weak magnetic fields 
strongly modify the convective stability criterion by causing heat to be conducted almost 
exclusively along magnetic field lines. This anisotropy in the thermal conductivity arises 
when the electron gyroradius is much less than the electron mean free path, a condition that 
is easily satisfied for realistic magnetic fields in most cases of interest. Balbus considered an 
equilibrium in which the magnetic field is in the xy- plane, and in which (3 — Sirp/B'^ ^ 1, 
where p is the pressure and B is the magnetic field strength. He showed that near marginal 
stability, the temperature of a rising fluid parcel is almost constant, essentially for two 
reasons. First, the parcel remains magnetically connected to material at its initial height. 
Second, near marginal stability a fluid parcel rises very slowly, so thermal conduction has 
enough time to approximately equalize the temperature along the perturbed magnetic field 
lines. As a result, the stability criterion becomes dT/dz > 0. When this criterion is satisfied, 
a rising fiuid parcel is cooler than its surroundings, and hence denser at the same pressure, so 
that it falls back down to its initial height. Parrish & Stone (2005, 2007) carried out numerical 
simulations that validated Balbus' analysis and extended it to the nonlinear regime. 

An immediate question arises, namely, why doesn't Balbus's stability criterion apply 
to stars? Although stellar plasmas are magnetized, heat is conducted in stellar interiors 
primarily by photons. As discussed by Balbus (2000, 2001), the conductivity is thus almost 
isotropic, so it is the Schwarzchild criterion that applies. The reason that the Schwarzchild 
criterion applies even if the isotropic conductivity is large is somewhat subtle. If ds/dz < 0, 
then adiabatic expansion would cause a slowly rising fluid parcel to be hotter and lighter than 
its surroundings, and hence buoyant. The effect of isotropic conductivity is then to relax 
the temperature in the parcel towards that of the immediately surrounding fiuid. However, 
because the conductivity is finite, the rising fiuid parcel's temperature is never decreased 
all the way to the temperature of its surroundings. The fiuid parcel thus remains slightly 
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hotter than its surroundings, and hence shghtly hghter at the same pressure, and the fluid 
is convectively unstable. Although isotropic conductivity does not modify the Schwarzchild 
stability criterion, it docs reduce the convective heat flux and the "efficiency of convection" 
in a convectively unstable fluid by decreasing the temperature difference 5T between rising 
fluid parcels and their surroundings (Cox & Guili 1968). 

Balbus's analysis has been extended in two ways by recent studies. First, Chandran & 
Dennis (2006), (hereafter referred to as [CD06]), investigated how the stability criterion is 
affected by the presence of cosmic rays that diffuse primarily along magnetic fleld lines. Like 
Balbus (2000, 2001), they assumed that 3 ^ 1 and took the equilibrium magnetic field to be 
in the a;j/-plane They showed that near marginal stability, the cosmic-ray pressure is nearly 
constant within a rising fluid element. This is because the fluid element remains connected 
to material at its initial height, and because fluid elements rise very slowly near marginal 
stability, so that there is plenty of time for cosmic-ray diffusion to approximately equalize the 
cosmic-ray pressure pcr along the perturbed magnetic fleld lines. [CD06] showed analytically 
that the stability criterion in the presence of cosmic rays is nkBdT / dz + dpa/dz > 0, where 
n is the total number density of thermal particles. 

More recently, Quataert (2007) considered buoyancy instabihties in a low-density, high- 
/3 plasma in the absence of cosmic rays, but allowing the equilibrium magnetic field to have 
a component in the z direction, parallel or antiparallel to the direction of gravity. Since the 
temperature is a function of 2;, the z component of the equilibrium magnetic field leads to 
an equilibrium heat flux. Quataert (2007) showed that this heat flux causes the plasma to 
become convectively unstable even if dT/dz > 0, so that the plasma is always convectively 
unstable if the magnetic fleld has a nonzero z component and a nonzero component in the xy 
plane. This heat-flux-buoyancy instability arises because of the geometry of the perturbed 
magnetic fleld lines in the plasma. For example, when the magnetic fleld is in the z direction 
and a fluid element is displaced upwards at a 45-degree angle with respect to the z axis, 
fleld lines converge as they enter the fluid element from "above" (i.e., from larger z). As a 
result, if dT/dz > then the parallel heat flux converges within the fluid element, causing 
the fluid element to become hotter than its surroundings, and thus less dense at the same 
pressure. Buoyancy forces then cause the upwardly displaced fluid element to rise unstably. 
(Quataert 2007) The nonlinear development of this instability was investigated numerically 
by Parrish & Quataert (2007). 

One of the open questions in this area of research is whether the buoyancy instabilities 
identifled in these previous studies for the /3 3> 1 regime still operate when the magnetic 
fleld strength is increased to the point that /3 < 1. We address this question in this paper. 
We consider the equihbrium geometry investigated by Balbus (2000, 2001) and [CD06], in 
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which the magnetic field is in the xy-plane - in particular, we set Bq = -Boy- We also allow 
for cosmic rays that diffuse along magnetic field lines, but now we allow f3 to take any value. 
We focus on wave vectors in the "quasi-interchange" limit, in which \kx\ ^ \ky\, {k^l ^ \kz\, 
and \kxH\ 3> 1, where H is the density scale height. This is the most unstable wave- vector 
regime for stratified adiabatic plasmas, because a small ky reduces the stabilizing effects of 
magnetic tension and a large k^ allows a rising fluid element to easily get out of the way of 
the next rising element beneath it by moving just a small distance in the x direction. (Parker 
1967, Shu 1974, Ferriere et al. 1999) We show analytically that the stability criterion in this 
limit is 

dT dp,. ^ dB' 

nks-r- + -1— + ^ J- > 0, (1) 

dz dz Sir dz 

and we present a heuristic derivation of this stability criterion from physical arguments. We 
also derive approximate analytical solutions to the dispersion relation for small-amplitude 
perturbations to the equilibrium in different parameter regimes, and compare these solutions 
to numerical solutions of the full dispersion relation. 

Our results are important for determining the convective stability of galaxy-cluster 
plasmas, in which cosmic-rays are often produced by central radio sources. Convection in 
intracluster plasmas is of interest because it may provide a mechanism for regulating the 
temperature profiles of galaxy-cluster plasmas and offsetting radiative cooling, thereby solv- 
ing the so-called "coohng-fiow problem." (Chandran 2004, 2005; Parrish & Stone 2005, 2007; 
Chandran & Rasera 2007). Our results are also important for determining the conditions 
under which the Parker instability can operate in the interstellar medium (ISM). Previous 
treatments of the Parker instability assume an adiabatic thermal plasma (Parker 1966, 1967, 
Shu 1974, Ryu 2003). Our results show that anisotropic thermal conductivity makes the ISM 
more unstable to the Parker instability, so that the instability can operate under a wider 
range of equilibrium profiles than was previously recognized. 

The remainder of this paper is organized as follows. In section [2] we outline the deriva- 
tion of the general form of the dispersion relation. In section [3] we specialize to the quasi- 
interchange limit, present our derivation of the necessary and sufficient condition for con- 
vective stability, and describe the properties of the unstable eigenmodes in plasmas that are 
very close to marginal stability. In section H] we present a heuristic, physical derivation of 
the stability criterion. We discuss the implications of our work for galaxy-cluster plasmas 
and the interstellar medium in sections [5] and O respectively. In section [7] we summarize 
our results, and in appendix |X] we present approximate analytic solutions and numerical 
solutions to the dispersion relation. 
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2. The general dispersion relation 



We begin with a standard set of two-fluid equations (Drury & Volk 1981, Jones & Kang 
1990), which we modify to include thermal conduction along the magnetic fleld: 



dp 

H 
dv 

Itt 

m 

dp 

'dt 

dpcr 
dt 



-pV ■ V, 

1 / 1 

-V U + Pcr + ^ +— B- VB + g, 

p \ ovr / Anp 



-BV v + B Vv, 

-7pV ■ V + (7 - 1) V 



bKii b ■ VT 



-7crPcrV ■ V + V 



hD\\ [h-Vpc 



(2) 
(3) 

(4) 
(5) 
(6) 



where d/dt = d/dt + \ ■ V, and where p is the plasma mass density, v is the velocity, p 
is the plasma pressure, Pcr is the cosmic-ray pressure, 7 is the ratio of specific heats of the 
plasma, 7cr is the effective ratio of specific heats for the cosmic rays, B is the magnetic 
field, b is a unit vector in the direction of the magnetic field, k\\ is the thermal conductivity 
along the direction of the magnetic field, D\\ is the cosmic-ray diffusivity along the direction 
of the magnetic field, and g is the gravitational acceleration. We have ignored cross-field 
conduction and diffusion in equations ([5]) and ([6]) since the gyroradii of the thermal particles 
are small compared to their Coulomb mean free path, and the gyroradii of the cosmic rays 
are small compared to the mean free path for cosmic-ray scattering. Equations ([2]) — ([6]) are 
closed via the equation of state for an ideal gas: 



p = Cv{i-l)pT. 



We take 



and consider an equilibrium in which 



Br 



-gz, 



Boy. 



(7) 



(8) 



(9) 



All equilibrium quantities (denoted with a "0" subscript) are taken to be functions of z only, 
and we set vq = 0. These assumptions lead to a condition for hydrostatic equilibrium of the 
form: 
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To introduce perturbations, we represent the variables in our two-fluid equations as sums of 
an equilibrium value and a small fluctuating quantity as follows: 

p = po + 5p (11) 
p = Po + 5p, (12) 



We employ a local analysis, in which we take the fluctuating quantities to be proportional 
to e*(^-^-^*), with 

kH:^l, (13) 



where 



H = 



d\n po 



dz 



(14) 



is the density scale height, which we take to be comparable to the length scales over which 
each of the equilibrium quantities varies. Substituting equations (fTTj) and (fT2|) . and analogous 
expressions for pcr, B, v, b, and T into equations (E]) — ([6]), and into equation (JTj), we obtain 
the following equations for the fluctuating quantities: 



. Sp ^ ^ dlnpo 

— too h «k • OV + OVz ; = U, 

Po dz 

• r Sp^ Sptot , 2 \d\nBo5B, 
tujov = —g — z — tk h Vj^ 

Po Po 



dz 5, 





5B d\nBo ...... 

I'Oi!—— = ikyO\ — OVz — J y (*k • ovj 

-DO 



y + iky— 

-Dn 



— lUJ 



5p 5p 
Po ^PoJ 



+ 5v, 



d In {pp 

dz 



■7 



dz 

-Ocond 



— luj h OVz ; \- ^7crk ■ ov = D\\ 



Pcr,0 



dz 



d\nTo6Bz ,2^'^ 

2 >-\Fcr 
>cr,0 



dlnpcrflSBz 

dz Bn 



,6pc 



6p 6p ^ 6T 
Po Po To ' 



where, 



and 



cond 



(7 - 1) i^\\To 
Po 



B^ 

Ptot = P + Pcr + 



(15) 

(16) 
(17) 
(18) 
(19) 
(20) 

(21) 
(22) 
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Equations (|T5|l - fl20|) may be reduced to an expression of the form: M ■ 5v 
a 3 X 3 matrix. For non-trivial solutions of this equation, we require |M| 
obtain the dispersion relation: 



Aquj^ + ^20;^ + A^u"^ + = 0, 



where, 



A2 
^4 

A, 



-e {u^ + vl) - klvl + g 



dhi Po 
dz 



kyvl {2u' + vl) - {kl + kl) 



+ {u^ + vl) g 



(ilnpo 

dz 



-k'klvW+{kl + k'^ (g' + u'g 



d\n Po 

dz 



0, where M is 
0, whence we 

(23) 



(24) 
(25) 

(26) 
(27) 



and 



u 



1 

Po L" 



Po 



UJ + 17] J UJ + IV 

and where in equation (1281) we have introduced the quantities 



T] 



V 



vl 



iTipO 



cond; 



(28) 



(29) 
(30) 

(31) 



where v^ is the Alfven speed, and t] and v are, respectively, the rates at which temperature 
fluctuations and cosmic-ray-pressure fluctuations are smoothed out along the magnetic field. 
Equations fl23|) - fl28|) represent the same result as that presented in equations (26) and (27) 
of [CD06] and we shall henceforth refer to this result as the "general dispersion relation." 
As we shall see, this relation constitutes an eighth-order polynomial equation in a = —iu, 
(where the change of variables is made so as to make all of the polynomial coefficients real). 
It is worthwhile to note that in the absence of cosmic rays and thermal conductivity. 



Po 



(32) 



where Cg is the adiabatic sound speed, and that if we take this limit, together with the 
limit of no stratification and g 0, the general dispersion relation reduces to the well- 
known dispersion relation obtained in ideal MHD. Thus the normal modes described by 
equation fl23|) may be viewed as modifications of the Alfven mode and the fast and slow 
magnetosonic modes of ideal MHD. 
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^3 



UJ4: = 



We now present the definitions of a number of frequencies that allow us to write the 
polynomial form of the dispersion relation more compactly. These are: 

= > (33) 

Po 

cul = klvl, (34) 

= f^g'sm^e, (35) 

2 ^sin^^ d , . 
^1 = ^^^^ 

= gsin^e-^lnT, (37) 
ln(p,,,oPo^"), (38) 

Tcr 

g'sin^6'^lnpcr,o, and (39) 

col = gsm'e^JnBl (40) 

(41) 

sin^^ = ^^. (42) 

The quantities u!\ and are the squares of the Alfven and isothermal sound-wave frequencies 
respectively. The quantity uf is the square of the usual Brunt- Vaisala frequency for buoyancy 
oscillations in the limit of vanishing cosmic-ray pressure and magnetic field. As we shall see 
below, the quantities c<j| and cj| serve to modify the frequency of these oscillations when 
the cosmic-ray pressure and magnetic field are non-vanishing. The quantities a;| and a;| are 
related to ujf and through the identities: 

u^l = iu;l + {j-l)gsin'9^, (43) 

2 2 , • 2/)^^"P0 /-,x 

<^4 = 7cr<^3 Tors' sm 9 . (44) 

2 



where we have defined 



We also define the quantities W and C: 

W' = ^i + XUjl + ^u;l (45) 

and, 

C^u;l + (46) 
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where in equation (H5|) . x = Pct,o/po- Finally, noting the identity 



— g sin^ 9 



d In po 

dz 



^0 - + 



we find that the general dispersion relation may be written 



a^a + aiO" + ■ ■ ■ + aycr 



where 



ao 
as 

04 



05 



06 

O7 



(z/ + r/) u,^ 



2 

s ' 



2 (i In Po 



^(7 + |)+.(l + x7c. + |l + 
1 2\ / 2 ^Inpo \ -2 



d\n Po 

dz 



([{1 - 1) + X7cr] + 1^ ^0 + (7 + X7cr) + (^7 + X7cr + |) 



7 + ^ ) (^0 + - ^0 + 7^1 



+ 



7/ 



/5 



1 + X7cr + -5 (^0 + C) - ^0 + (1 + X7 



+ [(7 + X7cr) (^0 + - ^0] 



^a{{^ (7 - 1) + XlcrV] ^0 + [z/7 + ^7 (1 + X7cr)] , and 



(47) 
(48) 

(49) 
(50) 

(51) 
(52) 



(53) 



(54) 

(55) 

(56) 
(57) 



We assume that \dlnpo/dz\, \d\npcr.o/dz\, and \dlnBQ/dz\ are of order H^^. We may 
thus conclude from equation (ITOi) that g ~ Ptot,oH^^ / po- We also assume that ptot is not 
much greater than p. Our assumption that \kH ^ 1| then allows us to write that 



(58) 



This inequality enables us to drop the u^'^gdliapo/dz terms in equations flSTj) . fl52|) . and 
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As a check on the results of this section, we show in appendix [B] that equation fHHl) 
reduces properly to the results obtained by Parker (1966, 1967) and Shu (1974) when cosmic- 
ray diffusivity is taken to be infinite and thermal conduction is negligible, and when the 
results of Parker (1966, 1967) and Shu (1974) are considered in the short-wavelength limit. 



\KH\ 


> 


1, 




\kx\ 


> 


\h 1 
\i^y\, 






> 




and 






1 





3. The quasi-interchange limit 

The most unstable modes in a gravitationally stratified adiabatic plasma threaded by a 
horizontal magnetic field are those for which \kx\ is very large, so that 

(59) 
(60) 
(61) 
(62) 

(Parker 1967, Shu 1974, Ferriere et al. 1999). We conjecture that the same is true when 
thermal conduction is taken into account, at least when the equilibrium magnetic field is 
horizontal, and thus we focus on this limit, which we call the "quasi-interchange limit." 

For very large \kx\, one set of modes consists of high-frequency magnetosonic-like waves. 
In the 3> 1 limit, these waves are stable [CD06], and we assume they are stable here as well. 
[We note, however, that in the presence of an equilibrium heat fiux (i.e. Bq^ 0), anisotropic 
conduction can cause magnetosonic waves to become overstable (Socrates, Parrish, & Stone 
2007).] To filter out these high-frequency waves, we assume that 

a -C ujg. (63) 

We also assume that \kx/ky\ is sufficiently large that 

< 1, and (64) 

< 1, (65) 



UJa 
V T] 



and that \kxH\ is sufficiently large that 



< 1; i = 0,...,5. (66) 



Using these inequalities and equation (1581) . we can rewrite the general dispersion relation as 
a 6th-degree polynomial equation, 

boa^ + bia^ + h^a^ + 63^^^ + h^a"^ + &5(^ + &6 = 0, (67) 
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where 



bo 
bi 
b2 
b. 



&5 

b. 



7 + XI cv + 



/3' 



M7 + |)+r/(l + X7c. + | 



(7 + Xlcr) {i^l + + C)-ujI+'^ {ujI + C) + ur](^ + 1^ 
[V (1 + Xlcr) + ^^7] K + ^0 + C) + (z/ + r^) 



(7 + X7cr) (^0 +C) -UJ, 



{ujI + C) [z/7 + r/ (1 + X7cr)] - (z/ + r/) 



and 



vrjuj\C. 



(68) 
(69) 
(70) 
(71) 
(72) 

(73) 
(74) 



3.1. Stability criterion 

To obtain the stabihty criterion for the modes described by equation fl67j) . we use the 
Routh-Hurwitz theorem, [see for example Levinson & Redheffer (1970)] To apply this theo- 
rem, we construct the matrix TZ from the (real) coefficients of the polynomial in equation ([6] 
where 

&5 
&4 
^3 
&2 

bi 
bo 



n 



(b, 
bo 




V 



&3 

&1 

bo 







bo 
b, 

b3 
62 







bo 

64 



\ 






bo 



(75) 



The Routh-Hurwitz theorem then states that for the real parts of the roots of equation ([6j 
to all take on negative values, it is a necessary and sufficient condition that the determinants 
of the principle minor matrices Aii of TZ all be positive-definite. This necessary and sufficient 
condition is the stability criterion for our plasma. The determinants of the principle minors 
of 71 are: 

det(l) = 61, (76) 
bi &3 

^0 &2 



det(2) 



(77) 



- 12 - 



det(3) 



det(4) 





61 63 65 






60 &2 ^4 






bi 63 




bi 63 fes 


bo &2 &4 ba 


61 63 6, 




60 62 &4 


&i 


&3 &5 





bo 


62 ^4 &6 








61 63 65 








60 &2 &4 


be 





61 63 


b. 



det(5) 



and 

det(6) = \n\ . 

After some algebra, wc find that these determinants may be expressed as 



det(l) 



/9 



2 



det(2) = 611/77 (^1 + -J + J (^a;o > ^ 
det(3) = 61 Ja;^ (0;^ + W"") + (^u;^ + + i/ryX^i (^a;^ + 



P 



det(4) = |j^a;iu;o^[iy^ + a;o^]V 



2 



(78) 



(79) 



(80) 



(81) 

(82) 
(83) 



(84) 



+ r,)ujl{ujl - ujlf + -K [{W^ + u^lfujl + Cc^oM + ^i(^i - ^0)'] 



det(5) 



+ {ur)ybiK 
2 



(l + |)c.o^C + |(c.o^-c.i)^ 



(85) 



2, .2 



212 



X 



r 2 

I (1/77)^ 61 + {ur]) JK biC + ujl{iy + r)) + k (o;^ + uol) + -ujI (1/ + 77) 

+ J^u;! [Xa;^ + (X + 1/ + 77) C] } and 



(86) 
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det(6) = |7^| = fee det (5) . (87) 
The quantities J and K appearing in the above expressions are defined as 

J= (7-l)r/ + x7cri^, (88) 

and 

ir=(7-l)z/ + x7crr/, (89) 

and are always positive. 

We first consider the case ky ^ 0. In this case, J, and the first two determinants are 
seen to be composed of sums of positive-definite quantities and so are themselves positive 
definite. By inspection, det(3) through det(6) are positive if C > 0, and thus C > is a 
sufficient condition for stability. On the other hand, equation (l87|l shows that if C < 0, then 
either det(5) or det(6) is negative. Therefore, C > is also a necessary condition for stability. 
If we fix the wavevector k, taking ky ^ the necessary and sufficient condition for modes 
at that k to be stable is then 

C > 0. (90) 

Since C = + kyv\, the smallest value of C is obtained in the limit ky — > 0. The necessary 
and sufficient condition for the plasma to be stable at all wavevectors in the quasi-interchange 
limit is thus 

> 0. (91) 

Using the definition of W"^ given in equation (H5|) . and the definitions of the frequencies ci;|, 
ul, and u;| given in equations (137|) . (!39l) . and (140|) respectively, we can rewrite equation (I9T1) 

as 

, dT dp^, 1 dB^ ^ 
dz dz Sir dz 



where we have dropped the zero subscripts on the equilibrium quantities. Equation ( l92l) 
shows that an "upwardly decreasing" temperature, cosmic-ray pressure, or magnetic pressure 
is destabilizing. 

We next consider the special case, ky = 0, which corresponds to pure interchanges. In 
this case, equation (p71) leads to the two non-trivial solutions, 



. = ±^-1, (93) 

where b is what remains of the coefficient 62 at A;^ = 0. By inspection we see that the 
necessary and sufficient condition for these modes to be stable is 6 > 0. For a vanishing 
cosmic-ray pressure, the condition 6 > reduces to 

- ^ > ^ (94) 
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where we have again dropped the zero subscripts on the equihbrium quantities. Equation flMl) 
is the result of Tserkovnikov (1960) for the pure interchanges as quoted by Newcomb (1961). 
The criterion > obtained above for the case fcy 7^ is more restrictive than the 
condition b > 0, since 7 > 1. Thus, W"^ > is the necessary and sufficient condition for the 
plasma to be stable to all modes in the quasi-interchange limit, including those with ky = 0. 



3.2. Eigenmodes near marginal stability 

In this section we consider the properties of unstable modes very near to the limit of 
marginal stability. We assume that ky 0, but take the limit kyH <^ 1 — that is, the parallel 
wave length is much longer than the scale height. Near marginal stability, the quantity feg 
in equation (!67|) approaches zero. There thus exists a solution to the dispersion relation 
in which a also approaches zero, for which the terms proportional to through in 
equation (1^ can be neglected. This solution satisfies the approximate equation 

a ~ (95) 

Making use of the fact that a for this mode, we can return to the results of section [2] 
(in particular, the equation M ■ (5v = 0) and show that to leading order in kyH 

kjv^ ik^p 

(96) 



and 



so that 



kySvy pg 

kzSvz ikzP 



ky5vy pg 



(97) 



\kx8vx + kz5vz\ <C \ky5vy\. (98) 



Thus, for modes with \kyH\ <^ 1 near marginal stability, most of the compression or expan- 
sion of the plasma occurs in the direction of the magnetic field rather than perpendicular to 
the magnetic field, despite the fact that ky is very small. We discuss the importance of this 
result further in the next section. 



4. Heuristic derivation of stability criterion 



In this section, we present a way of understanding the stability criterion in equation ([T]) 
in physical terms. We consider the same equilibrium discussed in section[2], in which g = —gi, 
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and Bo = -Boy, and we again take the plasma to be perfectly conducting, so that magnetic 
field lines are frozen-in to the fiuid. However, we now assume that the equilibrium is very 
close to marginal stability. We then imagine some mode in the plasma that causes a long 
and narrow magnetic fiux tube to rise upwards, as depicted in Figure [H For simplicity, we 
assume that the ends of the fiux tube are anchored at the fiux tube's initial height. We take 
the fiux tube to be very long, so that magnetic tension forces are very weak. Because the 
medium is arbitrarily close to marginal stability, the growth time or oscillation time for the 
mode is arbitrarily long. Thus, even though the fiux tube is long, there is plenty of time for 
conduction and diffusion to equalize T and pcr along the perturbed magnetic field lines. We 
assume that the total pressure, ptot = P + Pcr + B^/Stc, at each point along the flux tube is 
equal to the total pressure just outside the flux tube at that pointE 




We deflne An, AT, AB^ and Apcr, respectively, as the difference between the density, 
temperature, fleld-strength-squared, and cosmic-ray pressure at the highest point in our flux 
tube and the immediately surrounding medium, at a point in time when the top of the flux 
tube is a small distance Az above the flux tube's initial height. The constancy of T and Pcr 
along the flux tube yields the relations (accurate to flrst-order in Az/ H) 

AT = -Az^ (99) 
az 

and 

Apcr = (100) 

where dT/dz and dp^r/dz are the gradients of the equilibrium temperature and cosmic-ray 
pressure evaluated at the initial height of the flux tube. 



^Total-pressure variations are associated with high-frequency magnetosonic waves. These waves are stable 
at /3 » 1 [CD06], and we assume they are stable here as well. However, we note that Socrates, Parrish, & 
Stone (2007) have shown that magnetosonic waves can become unstable in the presence of an equilibrium 
heat flux, when Bqz 0. 
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Equations fl99l) and fllOOp tell us how to evaluate T and pcr in our flux tube. Evaluating 
in the flux tube is a little more involved. Assuming that the total pressure decreases with 
height, the fluid in the flux tube has to expand in order to achieve total-pressure balance. 
However, the manner in which the flux tube expands is not obvious. If the plasma expands 
primarily along the magnetic field, the cross sectional area of the flux tube will be constant 
along the flux tube, and thus so will the magnetic fleld strength. On the other hand, if the 
plasma expands perpendicular to the fleld, the magnetic fleld strength will decrease. Which 
type of expansion does the plasma favor? We answer this question analytically in section [221 
where we show that, near marginal stability, \ky6vy\ ^ {k^Svx + kz6vz\ for the low-frequency 
long-parallel- wavelength buoyancy instability in the large- |/Ca;.| limit. Thus, for this mode, 
most of the expansion V ■ v arises from the parallel motion. We note that this statement is 
stronger than the statement that \vy\ ^ \vx\, because we take \krc\ S> \ky\. 

How can we understand this result in physical terms? One way is by analogy to the 
SW analysis of the stability of ideal MHD plasmas, in the absence of thermal conduction. 
(Bernstein, Frieman, Kruskal, & Kulsrud 1958, Friedberg 1987) In this analysis, it is shown 
that if a mode expands in the direction perpendicular to the magnetic field, additional work 
must be done on the surrounding magnetic field. This requirement makes the mode more 
stable. To find the stability criterion, we must seek out the most unstable mode, which in 
this case is a mode that keeps the cross- sectional area of the fiux tube constant. 

Taking the cross-sectional area of the fiux tube to be constant, we can treat as 
constant along the fiux tube. This allows us to write that 

jr>2 

AB^ = -Az^. (101) 
az 

The condition that the total pressure inside the fiux tube matches the total pressure outside 
the fiux tube can be written as 



ksTAn + nksAT + Apcr + 



AB^ 



(102) 



0. 



Together, equations (!99|) through (11021) imply that 




(103) 



The stability criterion, equation ([T]), is thus the condition that the material inside an up- 
wardly displaced, long, and narrow fiux tube be denser than the surrounding medium. 
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5. Convection in galetxy cluster plasmais 

In many galaxy-cluster cores, the radiative cooling time is much shorter than the clus- 
ter's likely age (Fabian 1994). Nevertheless, high-spectral-resolution X-ray observations show 
that very little plasma actually cools to low temperatures. (Bohringer et al 2001; David et 
al 2001; Molendi & Pizzolato 2001; Peterson et al 2001, 2003; Tamura et al 2001; Blan- 
ton, Sarazin, & McNamara 2003). This finding, sometimes referred to as the "coohng-flow 
problem," strongly suggests that plasma heating approximately balances radiative cooling 
in cluster cores. 

A heating mechanism for cluster cores that has been studied extensively is heating 
by a central active galactic nucleus (AGN). The importance of such "AGN feedback" is 
suggested by the observation that almost all clusters with strongly cooling cores possess 
active central radio sources (Burns 1990; Ball, Burns, & Loken 1993; Eilek 2004) and by 
the correlation between the X-ray luminosity from within a cluster's cooling radius and the 
mechanical luminosity of a cluster's central AGN (Birzan et al 2004, Eilek 2004). One of 
the main unsolved problems regarding AGN feedback is to understand how AGN power is 
transferred to the diffuse ambient plasma. A number of mechanisms have been investigated, 
including Compton heating (Binney & Tabor 1995; Ciotti & Ostriker 1997, 2001; Ciotti, 
Ostriker, & Pellegrini 2004, Sazonov et al 2005), shocks (Tabor & Binney 1993, Binney & 
Tabor 1995), magnetohydrodynamic (MHD) wave-mediated plasma heating by cosmic rays 
(Bohringer & Morfill 1988; Rosner & Tucker 1989; Loewenstein, Zweibel, & Begelman 1991), 
and cosmic-ray bubbles produced by the central AGN (Churazov et al 2001, 2002; Reynolds 
2002; Briiggen 2003; Reynolds et al 2005), which can heat intracluster plasma by generating 
turbulence (Loewenstein & Fabian 1990, Churazov et al 2004, Cattaneo & Teyssier 2007) and 
sound waves (Fabian et al 2003; Ruszkowski, Briiggen, & Begelman 2004a,b) and by doing 
pdV work (Begelman 2001, 2002; Ruszkowski & Begelman 2002; Hoeft & Briiggen 2004). 

Another way in which central AGNs may heat the intracluster medium is by accelerat- 
ing cosmic rays that mix with the intracluster plasma and cause the intracluster medium to 
become convectively unstable. A steady-state, spherically symmetric, mixing-length model 
based on this idea was developed by Chandran (2004) and subsequently refined by Chandran 
(2005) and Chandran & Rasera (2007). In this model, a central supermassive black hole ac- 
cretes hot intracluster plasma at the Bondi rate, and converts a small fraction of the accreted 
rest-mass energy into cosmic rays that are accelerated by shocks within some distance Tgource 
of the center of the cluster. The resulting cosmic-ray pressure gradient leads to convection, 
which in turn heats the thermal plasma in the cluster core by advecting internal energy 
inwards and allowing the cosmic rays to do pdV work on the thermal plasma. The model 
also includes thermal conduction, cosmic-ray diffusion, and radiative cooling. By adjusting 
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a single parameter in the model (rgource), Chandran & Rasera (2007) were able to achieve a 
good match to the observed density and temperature profiles in a sample of eight clusters. 

The treatment of convective stability in the work of Chandran (2004, 2005) and Chan- 
dran & Rasera (2007) was based on the assumption that (3 = Sttp/B"^ ^ 1. The present paper 
investigates convective stability for arbitrary f3. One of the motivations for this work is the 
possibility that some clusters with short central cooling times ( "cooling-core clusters" ) may 
be in the /3 ~ 1 regime. For a fully ionized plasma with a hydrogen mass fraction X = 0.7 
and helium mass fraction Y = 0.29, 



Although many studies of the magnetic field strength in clusters of galaxies find B in the 
range of 1 — 5yuG (see, e.g., Kronberg 1994, Eilek & Owen 2002), some studies of Faraday 
rotation in cooling-core clusters find much stronger magnetic fields (Taylor & Perley 1993; 
Kronberg 1994; Taylor, Fabian, & Allen 2002). In the case of Hydra A, Taylor & Per- 
ley (1993) found a tangled magnetic field of ~ 30/iG, and Taylor, Fabian, & Allen (2002) 
found a tangled magnetic field of ~ 35/iG. The analysis of X-ray observations of Hydra A 
carried out by Kaastra et al (2004), when converted to a ACDM cosmology (see Chandran 
& Rasera 2007), indicate that Ue — 0.01 cm~^ and fc^T ~ 3.4 keV in Hydra A at r = 50 kpc. 
Equation fll04p thus shows that if B is indeed as large as 30/iG in the core of Hydra A, then 
j3 is of order unity. Values of /5 ~ 0.1 — 1 for cluster cores in several other galaxy clusters 
were reported by Eilek & Owen (2002). Although these studies suggest that /3 ~ 1 mag- 
netic fields could be common in cooling-core clusters, some caution is warranted here. Vogt 
& Ensslin (2005) have reanalyzed the Faraday-rotation data for Hydra A using an updated 
plasma density profile, and found an rms magnetic field of 7/iG, which corresponds to /5 = 15 
at r = 50 kpc in Hydra A. In the remainder of this section, we explore the implications of 
the condition /3 < 1 on convective instability in clusters, but the above uncertainty in the 
value of P in cooling-core clusters should be born in mind. 

In section [3] we showed that the necessary and sufficient condition for stability for a 
mode with fixed nonzero ky in the quasi- interchange limit {{k^l much larger than \ky\, \kz\, 
and H~^) is 



This equation shows that the magnetic field has two competing effects on convective stabil- 
ity. First, if the field strength decreases "upwards" (i.e., dB'^/dz < 0), the g(3~^ d In B'^/dz 
"magnetic-buoyancy term" in equation (11051) is destabilizing. On the other hand, the kyv\ 
"magnetic-tension term" is stabilizing. We can estimate the relative importance of the differ- 




(104) 




(105) 
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ent terms in equation fllOSp by defining the length scales Hf, Hb, and Hp via the equations 

d In Per 

(106) 
(107) 



dlnT ^ Per 



dz 



P 

d\nB 



dz 



and 



dz 

PI 
P ' 



The ratio of the magnetic-tension term to the magnetic-buoyancy term is then 



2gf3- 



klHBHp, 



(108) 



(109) 



while the ratio of the magnetic-tension term to the "fluid terms," g[d\nT/dz+{pcr/p)dlnpcr/dz] 
is 



(110) 



The magnetic field turns off the buoyancy instability at wavevectors for which the magnetic 
tension term dominates over both the magnetic buoyancy term and the fluid terms. At 
/? ~ 1, this happens for kyHpHs ^ 1 and kyHpHf ^1. If we take all the scale lengths 
to be comparable to the density scale height H, then at /3 ~ 1 magnetic tension turns off 



the instabihty for kyH S> 1, but is negligible for kyH -C 1. At (3 



1 and kyH 



1, the 



tension, buoyancy, and fluid terms are all comparable, and magnetic buoyancy and magnetic 
tension to some extent cancel out. When -C 1, magnetic tension dominates for kyH ^ 1, 
magnetic buoyancy dominates for kyH -C 1, and the two are comparable at kyH ~ 1, again 
assuming that Hp ~ Hb ~ Hf ~ H. 

To apply our results to galaxy-cluster plasmas, we imagine some hypothetical spherical 
equilibrium, and consider local modes at a radius r at some location where the radial compo- 
nent of the magnetic field vanishes, and where all the scale lengths are of order r. Our local 
analysis of a slab-symmetric equilibrium is strictly applicable only to modes with kyr ^ 1, 
that is, to modes with parallel wavelengths much less than the scale height. Our results 
show that such modes are stable when /3 < 1 because of the stabilizing effects of magnetic 
tension. 



6. The Parker instability in the interstellar medium 



The Parker instability is an unstable mode in a gravitationally stratified plasma that 
is driven by the buoyancy of the magnetic field and/or cosmic rays. (Parker 1966, 1967) 



- 20 - 



The Parker instability is thought to be important for the interstellar medium for several 
reasons. It has been argued that this mode, acting alone or in concert with the thermal 
instability (Field 1965), contributes to the formation of molecular clouds (Blitz & Shu 1980; 
Parker & Jokipii 2000; Kosinski & Hanasz 2005, 2006, 2007). It has also been suggested 
that the Parker instability is a mechanism for regulating the transport of magnetic fields and 
cosmic rays in the direction perpendicular to the galactic plane, and for driving the Galactic 
dynamo. (See, e.g., Parker 1992, Hanasz & Lesch 2000, Hanasz et al 2004). 

The Parker instability is very similar to the instability that we have investigated in 
this paper. Standard analyses of the Parker instability consider an equilibrium in which 
g = —gi,, B is in the a;|/-plane, p oc exp(— z/if), and if, T, /3, and Pcrlv are constant. 
(Parker 1966, 1967; Shu 1974, Ryu et al 2003). In early studies, the parallel cosmic-ray 
diffusion coefficient D\\ was taken to be infinite, since pcr was assumed to be constant along 
magnetic field lines. (Parker 1966, 1967; Shu 1974) On the other hand, Ryu et al (2003) 
considered the effects of finite Dy, as well as cosmic-ray diffusion perpendicular to magnetic 
field lines. All of these studies took the thermal plasma to be adiabatic. 

The analysis of the present paper extends our understanding of the Parker instability in 
two ways. First, we allow the equilibrium values of T, /3, and Pcr/p to vary with z. Second, 
we consider the effects of anisotropic thermal conduction. By doing so, we show that the 
condition dT/dz < makes a plasma more unstable to the Parker instability than when the 
equilibrium is isothermal. We also show that even if the equilibrium is isothermal, anisotropic 
thermal conduction makes a stratified plasma more unstable to the Parker instability than 
when the plasma is treated as adiabatic. The Parker stability criterion in the limit \krc\ oo 
for the equilibrium described above can be obtained from the ky ^ limit of equation (73) 
of Shu (1974): 

— + Per < il - 1)P. (Ill) 

ovr 

Multiplying this equation by —1/H and making use of the assumptions that dB'^/dz = 
—B^/H, dpcr/dz = —pcr/H, and dp/dz = —p/H, we can rewrite equation fillip as 

On the other hand, when anisotropic thermal conduction is taken into account, the stability 
criterion for this constant-temperature equilibrium from equation (JHSD is 

d fB^ \ , , 

+ Pcr >0. (113) 



dz \%Tl 

Since 7 > 1, equation (11131) is more restrictive than equation (11121) . and anisotropic thermal 
conduction allows for instability under a larger range of equilibria than when the plasma 
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is taken to be adiabatic. The reason for this is that as a fluid parcel rises and expands, 
anisotropic thermal conduction allows heat to flow up along the magnetic fleld lines into the 
rising fluid parcel. This heat flow increases the temperature of the rising fluid parcel relative 
to the adiabatic case and thereby lowers the density, making the fluid parcel more buoyant, 
as in the high-/? zero-pcr limit considered by Balbus (2000, 2001). 



In this paper we derive the stability criterion for local buoyancy instabilities in a strat- 
ifled plasma, with the equilibrium magnetic fleld in the y direction and gravity in the — z 
direction. We take into account cosmic-ray diffusion and thermal conduction along magnetic 
fleld lines and focus on the large-|/c2:| limit, which is the most unstable limit for adiabatic 
plasmas. Our work extends the earlier work of Balbus (2000, 2001) and [CD06] by allowing 
for arbitrarily strong magnetic flelds. Applying our work to galaxy-cluster plasmas, we flnd 
that increasing the magnetic fleld to the point that /? = Snp/B"^ < 1 would shut off buoyancy 
instabilities at wavelengths along the magnetic fleld that are much shorter than the equi- 
librium scale height. Our analysis also extends our understanding of the Parker instability 
by allowing the equilibrium values of T, j3, and Pcr/p to vary with z, and by accounting 
for anisotropic thermal conduction. We flnd that the interstellar medium is more unstable 
to the Parker instability than was predicted by earlier studies, which treated the thermal 
plasma as adiabatic. 

We thank Eliot Quataert for helpful discussions. This work was partially supported 
by NASA's Astrophysical Theory Program under grant NNG 05GH39G and by NSF under 
grant AST 05-49577 

A. Approximate analytical solutions and comparisons to numerical solutions 

In this appendix, in order to provide further insight into buoyancy instabilities in the 
equilibrium described in section [21 we derive a set of approximate analytical solutions to equa- 
tion (167]) . To check our analytical solutions, we compare them with numerical solutions of the 
general dispersion relation [equation fHH]) ]. We begin by deflning the quantities uj^^^ and kcrit'- 



7. 



Conclusion 




Po 



(Al) 



-22- 



and 



crit 



cond 



(A2) 



We assume that -Dcond ~ -Dy, and that the equihbrium cosmic-ray pressure is not very 
small compared to the equilibrium thermal pressure. The frequency uj^ei is comparable to 
the frequency of buoyancy oscillations in the medium except when the plasma is near the 
stabihty boundary for one of the buoyancy instabilities. The quantity fccrit is defined so 



that when L, 



we have Ur. 



vy — rucrit wc iiavc {ju^ei = 7] ^ v . These definitions and observations allow 
us to define two limits for which it is possible to derive approximate analytical solutions 
to equation fl^ : the "long-parallel- wavelength" hmit {\ky\ <^ fccrit) for which the diffusive 
frequencies v and r] are small compared to the buoyancy frequencies, and the "short-parallel- 
wavelength" limit {\ky \ ^ fccrit) for which the diffusive frequencies are large compared to the 
buoyancy frequencies. In the following subsections, we consider each of these limits in turn, 
and as a further check we also compare their high—/? limits with the results of [CD06]. 



A.0.1. Long-parallel-wavelength limit 

As stated above, in the long-parallel- wavelength limit we have rj,i> cu^gj, while all 
of the other quantities listed in equations (!33|) — (HO!) that remain in the dispersion relation 
[equation (l67|l ] are ~ u'^^f. We take advantage of this and apply the method of dominant 
balance (Bender & Orszag, 1978). We set a = cUrefC*", ujq = dJi-eft^o; = i^ref'^i, = ^ref<^2, 

UJ3 = UJj-e{UJ3, 004 = ^5 = ^vel^b, = ^rci^A, V = ^^relVi ^ = Cl^ref^', wherC € <^ 1, and 

substitute these into equation fl671) leading to the result: 



60^'^ + + 62^* + h^^ + 64^^ + ho + fee = 0, 



(A3) 



where the coefficients are now given by: 



hi 

k 



7 + X7cr + 

(7 + X7cr) [^l + ^0 + C) - '^o + I ('^o e^i>V (^1 + | 



e< (1 + X7cr) + z>7] i^i + ^'o + C) + (u + f,) 



p 



CjI + C 



~,2 



(A4) 
(A5) 
(A6) 
(AT) 



(7 + X7cr) [ul + C] -U'f 



ul + C + -l^ul+C 



(AJ 
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6r = eujl 



^0 + C) [n + fl{l + X7cr)] + , (A9) 

be = e^vfjuXC. (AlO) 



We may now solve for the two sets of approximate solutions to equation (1A3P by assuming 
that one set will satisfy a ~ tUref while the other satisfies a ~ ecUref- We present each of these 
solutions in turn in the following subsections. 

Adiabatic-buoyancy modes 

We refer to the set of solutions satisfying a ~ cUref as the "adiabatic-buoyancy modes" since 
these are of the order of the buoyancy frequencies except when near the limit of marginal 
stability. @ To obtain expressions for these modes we set, 

a = ao + eai + . . . , (All) 

where 5"o represents the lowest-order part of a and 5"i is the first-order correction. Substi- 
tuting this expansion into equation (lASp . collecting terms that are of like-order in e, and 
requiring sums of terms of like-order in e to vanish separately we find for the lowest-order 
terms: 

aa^ + + c = 0, (A12) 



where 



7 + X7cr + |, (A13) 



b = ( 7 + X7cr + I 1 ^ 



2 

(7 - 1) + X7cr + 



'^o + (7 + X7cr)cI;i, (AM) 
<! (7 + X7cr) C + [(7 - 1) + X7cr] ^0 \ ■ (A15) 



Restoring the dimensions we write the solutions with the notation: 



<^o,±± - ±\l ^ , (A16) 

where a, b, and c are the dimensional analogs of a, b, and c respectively, and where the 
left-most ± subscript on cr shall refer to the ± symbol inside the radical while the right-most 



^We note for clarity that in the limiting case of high-/? discussed in [CD06], these modes were referred to 
as the "adiabatic convective/buoyancy modes." 
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refers to the ± symbol outside the radical. We note that the solution cro,++ is unstable when 
c < 0. More explicitly, 



(7 + X7cr) C + [(7 - 1) + X7cr] ^0 < 0, 



(A17) 



which holds only when C is sufficiently negative. We define the buoyancy frequency N 
through the equation 

-1 



7^1 + X7cr^3 + 



7 + X7cr + 



P. 



which is related to W'^ through the identity 



2 , [(7 -l)+X7cr] ^,2 _ (7 + X7cr + 2 //3) 

7 + X7cr ° 7 + X7cr 



(A18) 



(A19) 



In a high— /3 plasma in the absence of cosmic rays, reduces to the Brunt- Vaisala fre- 
quency for buoyancy oscillations in a gravitationally-stratified medium. With A^ defined as 
in equation (lAlSp we may rewrite the condition in equation (1A17P as 



(A20) 



l + 2/[/5(7 + X7cr)] 

which, in the limit of high—/?, reduces to the corresponding result obtained previously for 
these modes in [CD06]. 



Quasi-isothermal buoyancy modes 

We refer to the set of solutions which satisfy a ~ euj-cf as the "quasi-isothermal buoyancy 
modes," since these modes are of the same order as the frequencies t] and u in the long- 
parallel- wavelength limit. To obtain the quasi- isothermal buoyancy modes, we now set o"o = 
0, so that a = edi + O(e^). Substituting into equations flASp — flAlOp and retaining only 
lowest-order terms we find: 

aal + 60-1 + c = 0, (A21) 

where now, 

a = (7 + X7cr)C + [(7-l) + X7cr]^o'. (A22) 
b = fC + Kul, (A23) 
c = z/r/C, (A24) 



and where the definition of K is given in equation (|89|) and we have introduced the positive 
definite quantity: 

/ = 7^^ + ^ (1 + X7cr) . (A25) 
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For notational convenience we have again reverted to the dimensional form of our solution. 
The discriminant of equation (IA2ip can be expressed in the form, 



^2 _ = { [7z/ -r]{l + X7cr)] C + [(7 - 1) z/ - xicvV] + (7 - 1) X7cr [C + ^^o)^ 

(A26) 

which is non-negative so that the solutions, 



-h ± - 4ac 

= ■ 



are both real. 



To examine the stability of these solutions we first note that when all three coefficients 
of equation flA21l) are non-zero, both roots will be negative if and only if the coefficients also 
all have the same sign. The necessary and sufficient condition for this is ac > 0. Inspection 
of equations flA22p — flA24p indicate that if c is positive then both a and b are also positive. 
If instead, c < 0, then we must have a < for both roots to have the same sign, in which 
case 6 < as well and the resulting roots are again both negative. Thus an unstable mode 
results only when c < < a, or explicitly: 

z/r^C < < (7 + X7cr) C + [(7 - 1) + X7cr] cul (A28) 
Once again making use of equation flAlQI) we may rewrite this result as 

an expression which again reduces in the high—/? limit to the result obtained for these modes 
in [CD06] where they are referred to as the "quasi-isothermal convective" modes. 



Note that when both inequalities in equation flA28P are satisfied, the quasi-isothermal 



mode is unstable, while if only the first inequality is satisfied, the quasi-isothermal mode is 
stable and the adiabatic buoyancy mode is unstable so that we are guaranteed an unstable 
mode whenever C < 0. Since these results must hold for arbitrarily small finite values of ky, 
we conclude that these inequalities are consistent with our stability criterion W"^ > 0. 



Comparison with numerical solutions 

As a final check of the solutions, we present numerical solutions to the general dispersion 
relation [equation (BHl) ]. using a suitably chosen set of parameters for comparison with the 
predictions of equation flA16l) . For the particular set of parameters chosen, the unstable 
mode is the quasi-isothermal buoyancy mode. Results for this mode are presented for the 
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unstable quasi-isothermal buoyancy mode in the top row of figure [2] for the cases /? = 1 and 
13 = 100. For the entire range of wave number shown, we have ky ^ fccrit? so we do not 
indicate the location of fccrit in the figure. The analytical solutions in these figures are seen 
to compare well with the numerical results, and the growth rates become negative where 
the quantity C passes from negative to positive values as indicated by the intersection of 
the vertical and horizontal dotted lines, demonstrating that the solutions honor the stability 
criterion C > 0. For the P = 1 case, the mode is unstable in the range < /Cy < 3.16 kpc^^, 
and the maximum growth rate occurs where ky = ky^^^^ ~ 2.1 kpc~^ and takes the value 
c"+,max — 0.114 Myr~^, while for the P = 100 case, the range of unstable wave numbers 
is < A;^ < 0.316 kpc~\ ky^^.^^ ^ 0.21 kpc"^ and (T+,max ^ 5.44 x 10"^ Myr"\ These 
results illustrate that a dynamically significant magnetic field results in a significantly larger 
maximum rate of growth for the instability as well as a much larger range of unstable wave 
numbers for these modes. 



A. 0.2. Short-parallel-wavelength limit 



We now take up the case ^ fccrit which we recall is the limit in which the rates of 
diffusion and conduction are large relative to the buoyancy frequencies. Once again we set 

cja = toircfCUA, as before; but we scale the diffusive frequencies according to z/ = e~^z>, and 
T] = e^^rj. and we define the scaled coefficients 60 = &i = biu^ef, &2 

64 = 64l^ref> ^5 = h^^ref^ ^6 = h^^rcC ^e obtaiu: 

with: 

2 



65(7 + 60 = 0, 



(A30) 



bo 

h 
h 

k 
k 

ba 



7 + X7cr + 



,-1 



[fl (1 + X7cr) + H[^l + (l^l + C)+{v + fl) 
(7 + X7cr) (ujI+C) -UJ^ 



~,2 







) - 


+ e ^z/?7 



Col 



-2 2 n 
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(A31) 
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(A33) 
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As in the previous section we obtain a 6th degree polynomial, and we anticipate that the 
solutions will split naturally into two sets, one in which a is of order the diffusive frequencies 
and which we will call the "diffusive modes," and one in which a is of order the buoyancy 
frequencies which we will call the "isothermal buoyancy modes." We consider each of these 
sets of solutions in the following subsections. 



Diffusive modes 



We expand a according to: 



e% + ... 



(A38) 



where now 5"_i is the lowest-order part of d and ao is the lowest-order correction to o-\. 
We then substitute the expansion given by equation ( ]A38p into equation ( ]A30[) and collect 
like terms in e. The lowest-order terms are of order e~^. Collecting these and requiring that 
their sum vanishes separately we obtain the quadratic: 



acr^ ]^ -|- ha_x + c = 0, 



where 



7 + X7cr + 



X7c, + 



2 

'^1 



(A39) 

(A40) 
(A41) 
(A42) 



and where we have once again reverted to the dimensional form of our solution. The dis- 
criminant is 



h - Aac 



+ AVT] (7 - 1) Xlc 



(A43) 



and is always positive so that cr_i is always real. Since & > 0, the mode is always damped. 
This result is just the generalization to the case of arbitrary 13 of the diffusive mode referred 
to in section 4.2 of [CD06]. 



Isothermal buoyancy modes 

For the isothermal buoyancy mode we set cr_i = in equation flA38l) and substitute into 
equation flA30l) again. The terms of order then yield a quadratic equation for ctq, whose 
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solution in dimensional form is 




-b ± a/62 _ 



(A44) 



where this time 



a 




(A45) 



b 



(A46) 



c 



ujIC. 



(A47) 



In the high-/3 limit, equation flA44p reduces to the results of [CD06]. 



Comparison with numerical solutions 



Once again we compare our approximate analytical solutions to the solutions obtained nu- 
merically from equation fHHl) for a set of parameters chosen to ensure that the solutions 
presented fall within the bounds of the short-parallel- wavelength limit. As before we present 
two examples of the unstable (isothermal buoyancy) mode; one for the case (3 = 1, and one 
for the case jS = 100. These are shown in the bottom row of figure [21 For each mode, the 
location of the critical wave number kait is marked with a vertical line. It is again seen that 
within the limits of validity of our analytical expressions, {\ky\ ^ kcrit), the approximate 
results and numerical results are in good agreement. Also as before, we mark the location 
where C passes from negative to positive values and observe that these solutions again honor 
the stability criterion C > 0. Finally we note that for both cases, the range of unstable wave- 
lengths is the same as for the long-parallel- wavelength limit while, in the (3 = 1 case, the 
maximum growth rate occurs for ky^^i^^ ~ 2.10 kpc^^ where we find o"+,+.max — 0.625 Myr~^, 
whereas for the (3 = 100 case we find ky^^^^x — 0.112 kpc"^ and cr+,+,rnax — 0.0325 Myr^^. 
Here again we observe the importance of the magnetic field in determining both the rate of 
growth and the range of unstable wave numbers. Comparing these results to those in section 
lA.O.lt we also note that high rates of diffusivity lead to growth rates several times larger 
than when these rates are low. 



B. Relation to Shu's (1974) analysis of the Parker instability 



In this appendix we show that the general dispersion relation given by equation fl48l) 
reduces properly to the results obtained by Parker (1966, 1967) and Shu (1974) when these 
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results are restricted to the short- wavelength approximation (kH ^ 1) assumed throughout 
this paper. In the limit considered by these authors, cosmic-ray difFusivity is taken to be 
infinite and thermal conduction vanishes, whence 1] —>■ 0, and u oo. Taking these limits 
in equation fHSl) and once again neglecting the terms involving gdlnpo/dz, our dispersion 
relation reduces to 

Coa^ + C2(T^ + 04^2 + C6 = 0, (Bl) 

where the coefficients are now 



Co 

C4 
C6 



7 + ^ ) (^0 + C) - ^0 + 7^1 



and 



(B2) 
(B3) 

(B4) 
(B5) 



We compare this result to the expression obtained by Shu (1974) in the limit of no rotation 
and no shear which is given by equation (53) of that paper. To do this we must take account 
of the differences in notation, the fact that Shu's expression is given in dimensionless form, 
and most importantly that given the global equilibrium assumed in Shu (1974), the scale 
height, H, in that paper is consistent with our definition of H given by equation (fT^ . It is 
also important to note that Shu's perturbations include a multiplicative "envelope" function 
which is an exponentially decreasing function of z above the origin. To account for this the 
quantity "fc" as defined in Shu (1974) must be set to zero. With all of these requirements 
accounted for we find that Shu's result may be expressed in our notation according to: 



0, 



(B6) 



where the coefficients are 



So 



ik 



PoH J 



00. 



2 2 



2 g kl 



^ + ^]H + C)-< + l<-^l^l2- 7^1 



(7 + 7t^i 



13 Hk"^ 



ikz 1 
l^kH 



ikz 1 
~kkH 



(B7) 
(B8) 

and (B9) 

(BIO) 



Because k^ < k, the terms in equations (1B9P and (IBlOp that involve the factor 1/kH must be 
neglected in the short-wavelength limit. Additionally, the term in equation flB8l) involving 
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kz can be seen to be small compared to uj^ as follows: 

ikzPo ikz ikz 1 



PoHuj^ k^H k kH 



<1, (Bll^ 



and similarly, the term in equation flB9p involving k^ can be seen to be small by comparison 
to 'yuj^ according to: 

2 = T7 72 ~ 72T72 ^ y^'-^) 

where the second-to-last relation uses the fact that g ~ pqH^^/ pq. Thus the terms by which 
the coefficients Sj of Shu (1974) differ from the coefficients Cj obtained from our general 
dispersion relation in the Parker /Shu limit and given by equations ( 1B2I) — ( iBll) . are just those 
that must be neglected in the short-wavelength limit to which we have restricted ourselves 
in this paper, so that in this limit the results we have obtained match those of Shu (1974). 

We can informally recover the necessary and sufficient condition for stability in this 
limit by applying the method of dominant balance (Bender & Orszag, 1978) which we have 
described in appendix |X1 to the dispersion relation given by equation fIBip with the quantity 
ojI used to scale the solutions, and with the quantities and C taken to be small compared 
to tUg . Thus to lowest-order we find that solutions, to equation fIBip that are of order uj^ 
satisfy: 

<s + C24s + l^l^l = 0. (B13) 
The discriminant of this expression can be expressed in positive definite form as 

K - ^ulf + ^'-ul K + 7-.^) + K)^ (I) ' > 0, (B14) 

and one can see by inspection of equation ( ]B3[) that C2 is also a positive-definite quantity. 
Thus the solutions, 



are always real and negative, and the modes of order cro,s ~ ^^s are always stable. 



The potentially unstable modes ctq^u can be found by seeking solutions to equation fIBll) 
that satisfy cr^/cjg ~ e where e -C 1, while continuing to assume that cUq and C are also small 
compared to ujI- One finds 



:y-^[(7-l)^o^ + 7C] 



ao,u = ±. /— [(7 - 1) ujI + 7C], (B16) 
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from which we conclude that these modes will be stable if and only if 

(7-l)cu2 + 7C>0, (B17) 

which is identical to the necessary and sufficient condition given by Shu (1974) when the 
latter is evaluated in the hmit kH » 1. 
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Fig. 2. — Unstable solutions of the dispersion relation for the long-parallel- wavelength limit (first row) and 
the short-parallel- wavelength limit (second row) for the cases j3 — 1 (left panels) and /3 = 100 (right panels). 
The solid curves correspond to the approximate analytical solutions given in the text, the dashed curves 
correspond to numerical solutions of equation (|48p. the intesection of the dotted lines indicate the location 
where C = 0, and in the bottom row, the dashed-dotted vertical line shows the location where ky = k„it. 
(for the long-parallel- wavelength limit in the top row, fccrit falls outside the range of wave numbers shown.) 



